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Abstract. In this paper we explore enumeration problems related to the num- 
ber of reachable configurations in a chip-firing game on a finite connected graph 
G. We define an auxiliary notion of debt-reachability and prove that the num- 
ber of debt-reachable configurations from an initial configuration with c chips 
on one vertex is a quasipolynomial in c. For the cycle graph C n , we apply these 
results to compute a near explicit formula for the number of debt-reachable 
configurations. We then derive polynomial asymptotic bounds for the number 
of debt-reachable and reachable configurations, and finally provide evidence 
for a quasipolynomiality conjecture regarding the number of reachable config- 
urations. 



1. Introduction 

On a finite connected graph G, we place a nonnegative integer amount of chips 
at each of the vertices. A vertex can fire if it contains at least as many chips as 
its degree, in which case it sends one chip to each adjacent vertex. Such a move 
is called a chip-firing move. Chip-firing on graphs is known to have applications 
to numerous related areas of mathematics, such as tropical geometry [3], statistical 
physics [5J, and combinatorics @]- 

In this paper, we consider the notions of reachability and debt-reachability and 
related enumerative problems. A configuration of chips C is reachable from an- 
other configuration G if there is a sequence of chip-firing moves that transforms 
configuration G to configuration C . If we ignore the restriction that a vertex 
must contain at least as many chips as its degree, we obtain the related concept of 
debt-reachability (see Section [2] for a precise definition). 

For a fixed graph G with n vertices, in Theorem |4.1| we show that the number 
of configurations debt-reachable from a configuration with c chips on one vertex is 
a quasipolynomial in c with degree n — 1 and period dividing k(G). In Theorems 
|4.9| and |6.1| we show that both the numbers of configurations debt-reachable and 
reachable from a configuration with c chips grow asymptotically as 



(n-l)!«(G) 

Finally, based on empirical evidence we conjecture that the number of configura- 
tions reachable from a configuration with c chips on one vertex is a quasipolynomial 



in c for sufficiently large values of c (Conjecture 7.1 ). 

In addition to this, we also prove various properties about the number of such 
configurations for the case where G is a cycle graph C n . For n = 3, in Theorems 
3.1| and |3.3| we compute an explicit expression for both the number of reachable 
configurations and debt-reachable configurations. For n > 3, in Theorem |5.5| we 
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prove that the number of configurations debt-reachable from the configuration with 
c chips on one vertex is equal to 

l/c + n-l\ 

+ f(c mod n) 

n\ n — 1 J 

where /(c mod n) is a constant which depends only on the remainder of c modulo 

n. 

The rest of this paper is organized as follows. In Section[2]we provide background 
and definitions. In Section [3j we examine chip-firing on C3 and prove Theorems 
|3 - 1 1 and |3.3| mentioned above. In Section [4j we prove quasipolynomiality and derive 
an asymptotic bound to the number of debt-reachable configurations. In Section [5] 
we apply the results of Section [4] to the class of cycle graphs C n . In Section [6j we 
prove that the number of reachable configurations is asymptotically equal to the 
number of debt-reachable configurations. Finally, in Section [7j we state Conjecture 
|7.I| and provide some progress towards a proof. 

2. Background and Definitions 

Unless otherwise stated, when we refer to a graph G, G is finite, undirected, and 
connected with n vertices. We discount the possibility of loops or multiedges; G 
will always be a simple graph. We will also implicitly assume the prior existence of 
some labelling of the vertices of G from 1 through n. In accordance with common 
notation, we will let C n refer to the cycle graph on n vertices, K n refer to the 
complete graph on n vertices, and P n refer to the path graph on n vertices. 

The degree of a vertex in a graph is equal to the number of other vertices it 
shares an edge with. We will denote the degree of vertex i by dj. 

Definition 2.1. The adjacency matrix A (alternatively A(G)) of a graph G is the 
n by n matrix which has a I in the ith row and jth column if there is an edge 
connecting vertices i and j and a otherwise. 

Definition 2.2. The degree matrix D (alternatively D(G)) of a graph G is the n 
by n diagonal matrix whose ith clement along the main diagonal is equal to di. 

Definition 2.3. The Laplacian matrix L (alternatively L(G)) of a graph G is the 
n by n matrix equal to D — A. 

It is easily shown that, if G is connected, L has rank n — 1 with a zero eigenvalue 
of (1, 1, . . . , 1)*. We will therefore often let L' denote any of the n — 1 by n — 1 
matrices formed by removing one row and column of L (usually, we will assume 
without loss of generality that this is the last row and column). We call L' the 
reduced Laplacian matrix of G. 

Definition 2.4. A spanning tree of G is an acyclic connected subgraph of G with 
n vertices. 

We will let k(G) count the number of spanning trees of G. A famous theorem of 
Kirchoff states that [7]: 

Theorem 2.5 (Matrix- Tree Theorem). For a graph G, 

k{G) = \detL'\ 
We proceed to define the notion of chip firing on a graph. 
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Definition 2.6. For a graph G, we call an n-tuple of nonnegative integers (ci, C2, C3, . 
a chip-firing configuration on G with q chips on vertex i. If we remove the nonncg- 
ativity constraint and allow the Ci to be arbitrary integers, we call such an n-tuple 
a general chip-firing configuration. If at least one of the c$ < 0, then we say the 
configuration is in debt, or is a debt chip-firing configuration. 

Definition 2.7. We say that vertex i can fire if c,; > di. In this case, vertex i loses 
di chips and every vertex adjacent to vertex i gains 1 chip. We call this operation a 
legal chip-firing move. If we remove the condition that we must have Ci > di, then 
we call the resulting operation a general chip- firing move, or just a chip-firing move. 
For the specific case where Cj < , we will call the operation a de&i chip-firing move. 
A debt chip-firing move always results in a debt configuration. 

Note that the total number of chips on a graph stays constant under a chip- firing 
move. Moreover, as an operator on the space of chip-firing configurations of a graph 
G, firing vertex i acts by translating the current configuration by the zth row of 
L(G). That is, firing vertex i sends configuration C to the configuration C + Lei, 
where is the ith unit vector. A corollary of this is that chip-firing moves are 
abelian; firing vertex i and then firing vertex j will result in the same configuration 
as firing vertex j and then firing vertex i. Similarly, we also have the following 
consequence: 

Theorem 2.8. Firing each of the vertices once leaves the configuration invariant. 

Proof. This follows from the fact that (1, 1, . . . , 1)* is a zero eigenvalue of the Lapla- 
cian, as noted above. □ 

Definition 2.9. We say that a chip configuration C' is reachable from a chip 
configuration C if there exist a sequence of legal chip-firing moves which transform 
C into C' . We say that a chip configuration C' is debt-reachable from a chip 
configuration C if there exists a sequence of chip-firing moves (possibly debt chip- 
firing moves) which transforms C into C' . Note that both C and C' must be 
non-debt chip configurations. 

As a direct consequence of Definition |2 . 9| and the preceding observation, we have 
the following useful check for debt-reachability: 

Theorem 2.10. A configuration C is debt-reachable from the configuration C' iff 
there exists a solution x € Z n to Lx = C — C' . 

Theorem |2. 10| allows us to define the following: 

Definition 2.11. A vector x £ Z n is a firing vector for the transition from C' to C 
if Lx = C — C' . If x = (xi, %2, ■ ■ ■ j x n ), then Xi corresponds to vertex i being fired 
Xi times. Since the only zero eigenvector of L is (1, 1, . . . , 1)*, if Lx = C — C' has 
a solution, there is a unique solution such that min(ir,) = 0. We call this unique 
solution the reduced firing vector for the transition from C' to C . 

Finally, we define some terminology related to the theory of Ehrhart polynomials. 
Let an integral polytope be a polytope whose vertices have integer coordinates. 

Definition 2.12. The Ehrhart polynomial of an integral polytope P in R™ is de- 
fined to be: 



i(P,q) = #(qPnZ n ) 
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That is, i(P, q) is equal to the number of lattice points inside the q-fold expansion 
of the polytope P. It is a theorem of Ehrhart pQ that i(P, q) is indeed a polynomial 
of degree n in q for integer values of q. 

Likewise, let a rational polytope be a polytope whose vertices have rational 
coordinates. 

Definition 2.13. A quasipolynomial q(k) is a function from Z to Z of the form: 

q{k) = C d {k)k d + C d _ 1 (fc)fc d - 1 + • ■ • + C (fc) 

where the Ci{k) are periodic functions of k with integer period. We call the 
period of q the lowest common multiple of the periods of the Cj(fc). 

Then we also have the following definition. 

Definition 2.14. The Ehrhart quasipolynomial of an rational polytope P in W 1 is 
defined to be: 



i(P,q) = #(qPnZ n ) 



As in Definition 2.12 it is a theorem of Ehrhart [7] that i(P,q) is indeed a 
quasipolynomial of degree n in q for integer values of q. Moreover, if kP is an 
integer polytope, then the period of i(P,q) divides k. 

3. Enumeration of reachable configurations on C 3 

In this section, we explicitly compute the number of debt-reachable and reach- 
able configurations from the configuration (c, 0, 0) on the triangle graph C3. The 
techniques we use to do this are also illustrative of the techniques we will apply 
later in the paper to more general graphs. 

Theorem 3.1. The number of configurations that are debt-reachable from (c, 0, 0) 
is equal to: 



3k 2 - 


f 3fc- 


-2 




2 




3k 2 - 


f 5k- 


-2 




2 




3k 2 - 


f 7k- 




2 



if c = 3k, 



if c = 3fc + 1, and 



ifc = 3k + 2 

Proof. We wish to compute how many configurations are debt-reachable from (c, 0,0). 
Since we are only concerned with debt-reachability, we do not care about the or- 
der of the firings, so we can just assume that the firing vector for this transition 
is equal to (x\, X2, £3), for xi,X2,x^ G Z. After such an operation, we are in the 
configuration given by: 



(c — 2x 1 + x 2 + x 3 ,Xi - 2x 2 + x 3 , x 2 — 2a; 3 + X\) 
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In order for this configuration to be valid, all components must be greater than 
or equal to 0. This gives us the following set of inequalities. 



— 2xi + x 2 + x 3 + c > 
X\ — 2x 2 + x 3 > 

X'2 — 2x$ + Xi > 

Now, since firing vectors are equivalent modulo translation by (1,1,1) (from 



Theorem 2.8), set X\ = 0. This gives us the following set of inequalities in two 
variables: 



X2 + X3 > -c 
X3 — 2x2 > 
x 2 — 2x3 > 

These inequalities define a triangular region of the X2 — £3 plane with vertices 
at (0, 0), (— ^, — § ), and (— |, We wish to count the number of lattice points 

within this rational triangular region. The theory of Ehrhart polynomials tells us 
this has to be a quasi-polynomial of degree 2. To actually compute it, first note 
that if n is divisible by 3, we can use Pick's formula. If c = 3k, then the vertices 
of this triangle are (—2k, —k), (—k, —2k), and (0, 0). A simple count shows us that 
there are 3k lattice points on the boundary of this triangle. Moreover, the area of 
this triangle is equal to Since Pick's formula states that: 

/='+!-' 

where A is the area of the region, / is the number of lattice points in the interior, 
and B is the number of lattice points on the boundary. Applying this to our triangle, 
we have that: 

3fc 2 3k 

-T = I + Y- 1 

from which we have 

_ 3k 2 - 3k + 2 
~ 2 

The total number of lattice points in the triangle is thus: 

. D 3k 2 + 3k + 2 
I + B= 

Now, if c = 3k + 1, let us count how many lattice points are added to the 
triangular region (i.e., that lie in the triangular region for n = 3k + 1 but don't for 
c = 3fc). These are simply the lattice points which lie on the line X2+X3 — —3k— 1 
for — I > X2 > — Since X2 must be an integer, this implies that — (fc + 1) > x 2 > 
—2k. Since for each such X2 we can find a unique integer value of X3, we therefore 
have that for the case where c = 3k + 1, the total number of lattice points equals: 

3k 2 + 5fc + 2 
2 
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Likewise, when c = 3fc + 2, the number of additional points is simply the number 
of integral x 2 that satisfy — | > x 2 > — 4f , or equivalently, — (k + 1) > x 2 > 
— (2k + 1). Hence, there are k + 1 additional points, so the total number of lattice 
points for c = 3/j + 2 is: 

3fc 2 + 7fc + 4 
2 

Each of these corresponds to a unique debt-reachable configuration, so we are 
done. □ 

To compute the number of reachable configurations, we need the following result: 

Theorem 3.2. // a chip configuration (ci, 02,03) of C3 is debt-reachable from a 
configuration but not reachable from the same configuration, then (01,02,03) is a 
permutation of (c, 0, 0) (where c = Ci + c 2 + c 3/ ). 

Proof. It is clear that the configurations (c, 0, 0) (and its cyclic permutations) arc 
not reachable from any other configuration, because there is no possible prior con- 
figuration (any such prior configuration would either have the form (c + 2, —1, —1) 
or (c — 1,2,-1)). Assume that we are trying to reach some configuration C that 
is debt-reachable from configuration C. Let x — (x\, x 2 , £3) be the reduced firing 
vector of the transition from C to C . If X\ = x 2 = x 3 = 0, C and C are the 
same configuration, and we are done. Hence assume that at least one of the Xi is 
greater than 0, or alternatively that x\ + x 2 + x 3 > 0. We will induct on the value 
of x\ + x 2 + £3; that is, we will show that unless C — (c, 0, 0) (or a cyclic permu- 
tation) , we can fire a vertex of C and get to a configuration where the components 
of the reduced firing vector sum to x\ + x 2 + X3 — 1. 

As long as c > 3 (and for c = 1,2,3 we can check this lemma easily by hand), 
there will be one vertex with at least 2 chips on it. Without loss of generality, let 
this vertex be vertex 1. If x\ > 0, fire vertex 1, so we arrive at a new configuration 
C" where we have to fire vertex 1 x\ — 1 times, vertex 2 x 2 times, and vertex 3 £3 
times, and we have reduced this to a smaller case. Otherwise, we must have x\ = 0. 
Similarly, if either vertex 2 or 3 has at least 2 chips on it and if either x 2 or x 3 is 
greater than 0, fire the corresponding vertex and reduce the sum x\ + x 2 + x 3 by 
one, as before. 

Furthermore, if x 2 — 0, then only x% is greater than 0. We then claim that C3 
must have at least 2 chips on it, because otherwise the final configuration would 
be (ci + x 3 ,c 2 + x 3 ,c 3 — 2x 3 ), which is not valid if C3 < 2x 3 . Hence both C3 > 2 
and x 3 > 0, so we can fire vertex 3 and decrease X\ + x 2 + x 3 by one. Likewise, by 
symmetry we cannot have X3 = and x 2 > 0, so both x 2 and X3 must be greater 
than 0, and consequently, both c 2 and C3 must be less than 2. 

There are then at most four possibilities we need to consider; (02,03) = (0,0), 
(0,1), (1,0), and (1,1), or equivalently, (ci,c 2 ,c 3 ) = (c,0,0), (c-1,0,1), (c-1,1,0), 
and (c— 2, 1, 1). We can check that the only case where there are values of x 2 and x 3 
larger than which result in a legal chip-firing configuration is the case (c — 2, 1, 1), 
where if we fire both vertices 2 and 3, we end up at (c, 0, 0). This means, however, 
that C = (c, 0, 0), and hence we are done. □ 

This immediately leads to the following theorem. 
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Theorem 3.3. The number of reachable configurations from the inital configuration 
(c, 0, 0) on C3 is egual to the number of debt-reachable configurations when c is not 
divisible by 3, and is 2 less than the number of debt-reachable configurations when 
c is divisible by 3. Explicitly, the number of reachable configurations equals: 

3fc 2 +3fc-2 

— if c = 6k 



2 

3fc 2 + 5k + 2 
2 

3fc 2 + 7k + 4 



if c = 3k + 1 



if c = 3k 



Proof. We can check (via Theorem 2.101 that (0,c, 0) (and similarly, (0,0, c)) are 



only debt-reachable from (c, 0,0) when 3\c. Therefore the number of reachable 
configurations is just the same as the number of debt-reachable configurations when 
c is not divisible by 3, and is 2 less than the number of debt-reachable configurations 
when c is divisible by 3. □ 

4. DEBT-REACHABLE CONFIGURATIONS ON ARBITRARY GRAPHS 

In this section, we derive asymptotics and properties of the number of configu- 
rations debt-reachable from the configuration (c, 0, 0, . . . , 0) on an arbitrary graph 
G. 

Theorem 4.1. The number of configurations debt-reachable from (c, 0, 0, . . . , 0) on 
an arbitrary graph G is a quasipolynomial in c of degree n — 1 with period dividing 
k(G). 



Proof. Recall that (Theorem |2. 10 ) a configuration C' — (ci, C2, . . . , c„) is debt- 



reachable from C — (c, 0, 0, .., 0) iff C' — C = Lx for some vector x = (xi,x%, ■ ■ ■ , x n ) 
in Z n . Since (1, 1, . . . , 1) is the only zero eigenvector of L, we can impose the 
additional condition x n = 0. There is then a bijective correspondence between 
vectors x £ Z™ with x n = where C + Lx has nonnegative components and debt- 
reachable configurations C' . We therefore want to count the number of such x 
such that C + Lx has nonnegative components. Such x are determined by n linear 
inequalities in the n — 1 variables x\ through x n -\. Each of these inequalities 
determine half-spaces of K" _1 . Of these half-spaces, n— 1 pass through the origin. 
The remaining half-space is of the form a • x > — c for some vector a £ Z n . These n 
half-spaces thus determine a rational (n — l)-simplex in (n — l)-dimensional space 
with one vertex at the origin and n—1 other vertices at points cvi for some Vj £ Q" 
independent of c, for 1 < i < n — 1. 

Therefore, if we let P be the simplex formed when c = 1, the simplex formed 
for arbitrary c is just the c-fold dilate cP of P. Since the number of lattice points 
contained in cP is given by the Ehrhart quasipolynomial of P, the number of con- 
figurations debt-reachable from (c, 0, 0, . . . , 0) is a quasipolynomial in c. Moreover, 
since P is (n — l)-dimensional, this quasipolynomial has degree 

It remains to show that the period of this quasipolynomial divides n(G). To do 
this, it suffices to show that K,(G)vi £ Z' i_1 for each v*. Let L\ denote the reduced 
Laplacian matrix formed by removing the last column and the ith row from L. 
From the construction of v,-, we know that v,- satisfies: 
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L<v < = (l,0,...,0) t 



so 



v < = (£<)- 1 (1,0,...,0) < 
Since is an integer matrix, an elementary theorem of matrix algebra says 
that (det L^)(L9 _1 i s an integer matrix (namely, the cofactor matrix). Therefore 



(det I/j)vj is an integer vector. By Theorem 2.5 |deti^| = k(G), we have that 



k(G)Vj is an integer vector, as desired. □ 

Before we continue further, note that debt-reachability is an equivalence relation 
on the set of chip-firing configurations on a graph G with a total of c chips, since if 
configuration C is debt-reachable from configuration C , C is also debt-reachable 



from configuration C (this is an immediate consequence of Theorem 2.101. This 
equivalence relation therefore partitions this set into several subsets. We will call 
these subsets blocks, and let b c count the number of such blocks for the set of 
configurations with a total of c chips. 

We will now establish some properties about the sequence b c . 

Theorem 4.2. For all c > 1, b c > b c -\. 

Proof. Note that if configurations C — (ci, c 2 , . . . , c„) and C — (c[, c' 2 , ■ ■ ■ , c' n ) 
are not debt-reachable from each other (i.e. lie in different blocks), then the con- 
figurations D = [c\ + 1, C2, . . . , c„) and D' = (c[ + l,c' 2 , . . . , c' n ) are likewise not 
debt-reachable from each other. Also, if C and C are configurations with c — 1 
chips, then D and D' are configurations with c chips. It therefore follows that 
b c > bc-ij as desired. □ 

Theorem 4.3. For all c, b c < n(G), and for some sufficiently large C, be — k(G). 

Proof. Consider the integer lattice in n-dimensional space containing the points Lx 
for all vectors x in Z n (i.e. it is generated by the rows of L). We refer to this lattice 
as the chip-firing lattice, since an element (r 1; r 2 , . . . , r n ) lies in this lattice iff there 
is some sequence of chip-firing moves which sends a configuration (c\, . . . , c n ) to 
(ci + n, C2 + r2, . . . , On + r n ). We denote this lattice by CFL. 

Note that for any element r = (n, r<i, . . . , r n ) in this lattice, r% + r 2 + ■ ■ • + r n — 
0. We can therefore project this lattice onto the subspace where r n = and 
uniquely identify any element in the original CFL by the first n — 1 components, 
(ri, r2, . . . ,r n _i). We denote this modified lattice by CFL'. Furthermore, note 
that since the last row of L is an integer linear combination of the first n — 1 rows 
of L (since the sum of all n rows is zero), the lattice CFL' consists of all points 
i'x, where x lies in Z™ -1 and L' is the reduced Laplacian matrix. 

We then have that the number of elements of Z™ -1 /CFL' = Z" _1 /L'x is equal 
to det L' (this is a standard fact in the theory of lattices). In other words, there are 
at most det L' distinct "residue classes" of chip-firing configurations; each possible 
chip-firing configuration is debt-reachable from one of a particularly chosen set of 



detL' configurations. By Theorem 2.5 det L' = n(G), so b c < k(G) for all c 

Now, b c is actually counting the number of distinct residue classes of the lattice 
points in Z" _1 contained in the rational (n— I)-simplex defined by the n inequalities 
ri > — c, 7*j > for 2 < i < n — 1, and (r x + r 2 + ■ ■ ■ + r n _i) < 0. 
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This region scales with c, so eventually it will contain a sphere large enough 
that must contain a fundamental region of CFL' . At this point, it will contain a 
point for every single residue class, so there is some C such that be = k(G), as 
desired. □ 

Remark 4.4. The finite abelian group Z™ -1 /CFL' of equivalence classes modulo 
the lattice CFL' is commonly known as the sandpile group of the graph G. Much 
research into chip-firing on graphs has focused on the properties of these groups; 
see [H[S] for details. 



As a corollary of Theorems 4.2 and 4.3 we have 



Corollary 4.5. There exists a C such that b c = n(G) for all c > C . 



As another consequence of Theorems |4.2| and |4.3[ we have the following inter- 
esting (albeit slightly tangential) result: 

Corollary 4.6. If G is a tree, then all chip-firing configurations on G with c chips 
are reachable from eachother. 

Finally, we cite the following well-known lemma regarding the approximation of 
volumes. 

Lemma 4.7. Let P is a polytope in M n with content V(P), and let L be a lattice 
in W 1 whose fundamental domain has content n(L). Then as c tends to infinity, 
the number of points of L inside cP is asympotically equal to 

c n V(P) 

W) 

Proof. Refer to [J. □ 
We are now ready to derive an asymptotic estimate for the sizes of the blocks in 

b c . 

Theorem 4.8. The size of each of the blocks in b c is asymptotically equivalent to 

1 /c + n-l x 



k(G) V n-1 

as c approaches infinity. This in turn is asymptotically equivalent to 



(n- 1)!k(G) 



Proof. Note that there are a total of ( c ^™7 1 ) chip-firing configurations with a total 
of c chips on a graph with n vertices. Since there arc a total of k(G) blocks for 
sufficiently large c (by Corollary |4.5[ ) , it suffices to show that all blocks of b c have 
asymptotically the same size. 

Let P c be the polytope determined by the convex hull of the region defined in 
the proof of Theorem |4.3| Let B be a specific block of b c , and let v be a chip- 
firing configuration belonging to B. Then if we consider the (uncentered) lattice 
containing the elements v + L'x for x e Z" _1 , \B\ is simply equal to the number 
of points belonging to this lattice that also lie within P c . 

Now, as mentioned in the proof of ■ 



Hence, by Lemma 4.7 we have that 



4.3 the region P c scales with c (so P c — cP\). 



B\ asymptotically approaches: 
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c n - l V(Px) 
H{v + L'x) 

Since the volume of the fundamental region of a lattice is unchanged by translat- 
ing the lattice, //(v + L'x) is independent of v (in fact, it is equal to dct V = n(G)). 
Hence \B\ asymptotically approaches the same value regardless of the choice of block 
B, and therefore all blocks of b c have asymptotically the same size, as desired. □ 

But now, since (by definition) the configurations debt-reachable from a chip- 
firing configuration are those that belong to the same block as the configuration, 
we have the following main theorem as an immediate consequence of the preceding 
result. 

Theorem 4.9. The number of configurations debt-reachable from an arbitrary con- 
figuration with c chips asymptotically approaches 



(n- l)!/s(GQ 

as c approaches infinity. 



In particular, since by Theorem |4.1| the number of debt-reachable configurations 
from (n, 0, . . . , 0) is given by a quasipolynomial of degree n — 1, this implies that 

the leading term of this quasipolynomial is always equal to -, r-j — -r— . 

(n — ±)Ik(G) 

5. Debt-reachable configurations on C n 

In this section we apply the results of the previous section to obtain an exact for- 
mula for the number of debt-reachable configurations on C n for an infinite number 
of configurations of the form (c, 0, . . . , 0). 



First, note that as a corollary of Theorem 4.9 (along with the fact that fc(C n ) = 
n), we know that the number of debt-reachable configurations from (c, 0, . . . , 0) on 
C n is asymptotically equal to c n , . This matches up with our explicit formulae for 
C3, where the number of debt-reachable configurations is asymptotically equivalent 

To strengthen this result, we must make some additional observations regarding 
the numbers b c in the case of C n . Let represent the configuration where there 
is one chip on vertex i and zero chips elsewhere, and let ce^ similarly represent the 
configuration where there are c chips on vertex i and zero chips elsewhere. 

Theorem 5.1. For the case of G = C n , b\ = n. 

Proof. We must show that the configuration e, is not debt-reachable from ej for 
all i 7^ j. Due to the symmetry of C n , it suffices to just show that ej is not 
debt-reachable from e\ for i 7^ 1. We proceed to do this. 

Assume is debt-reachable from ei, and let x be a corresponding firing vector. 
We thus have that: 

Lx = e.j — e\ 

Writing out x = (xi,X2, ■ ■ ■ ,x n ), this is equivalent to the following set of n 
equalities: 
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x n - 2xi + x 2 = -1 
Xi-i - 2x t + x l+ i = 1 
Xj—i — 2xj + Xj + i = for j =/= 1, i 

Define y^ = x$+x — Xi, where indices are taken modulo n (so x n+ \ = xi). This 
system is then equivalent to: 

Vl-Vn = -1 

Vi - Vi-i = i 
yj-Vj-i = Oforj^l,i 
Note that we also have the condition y% + y% + ■ ■ ■ + y n = 0. This system is 
easy to solve (seeing as the yj take on only two different values), and has solution 
yj = — f for f < j < i — 1 and y 3 = otherwise. Since there is no value of 
i between 2 and n inclusive for which these yi are integers, as a consequence, the 
Xi cannot possibly all be integral, and therefore there is no such firing vector, as 
desired. 

□ 

Corollary 5.2. For G = C n , b c = n for all c. 



Proof. Follows from the fact that n(C n ) — n along with Theorems 5.1 and 4.3 □ 

Theorem 5.3. For C n , there exist no pair of distinct indices i,j such that cei is 
not debt reachable from cej iff c is relatively prime to n. 

Proof. We proceed identically to the proof Theorem |5.1| Such a pair of indices can 
only exist if there is an integer solution to: 



2/i - Vn = -c 
Vi ~ Vi-i = c 
Vj-Vj-i = Oforj^ 1 ^ 
with yi + 2/2 + • • • + Vn = 0. Similarly to before, the solution to this system of 
equations is yj — c( i ^i — 1) for 1 < j < i — 1 and yj = c( lz ^) otherwise. If c is 
relatively prime to n, there can exist no integer solution for yj by the same logic as 
On the other hand, if gcd(c, n) = r > 1, then if we set i = 1 + 



5.1 



in Theorem 

we have an integral solution for the j/j , as desired. □ 

We can now prove the main theorem of this section. 

Theorem 5.4. If c is relatively prime to n, then the number of configurations 
debt-reachable from the configuration (c, 0, 0, . . . , 0) on C n is equal to: 

1 /c + n - 1 

n \ n — 1 

Proof. Clearly, the number of configurations debt-reachable from the configuration 
cei must equal the number of configurations debt-reachable from the configuration 
ce,-, by the symmetry of C n . Theorem |5.3| shows that if c is relatively prime to n, 
all the configurations cej lie in distinct blocks, and theorem |5 . 1| show that there are 
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always exactly n blocks. Hence, since there are a total of (^"j/) configurations, 
each block must contain — ( c i™7 ) configurations, from which the result follows. 

□ 

We can use Theorem |5.4| to deduce the following information about the number 
of chip-firing configurations for values of c not relatively prime to n. 

Theorem 5.5. The number of configurations debt-reachable from (c, 0, 0, 0) on 
G n (for c not necessarily relatively prime to n) is equal to 



1 fc + n- I s 
n \ n — 1 

where f(c mod n) is a constant which depends only on the remainder of c modulo 



/(c mod n) 



Proof. By Theorem |4.1[ we know that this number of debt-reachable configura- 
tions is a quasipolynomial in c whose period we can take to be k(C„) = n. We also 
know that the number of debt-reachable configurations must be weakly increasing 
with c (adding a chip cannot possibly decrease the number of debt-reachable con- 



figurations). These two observations along with Theorem 5.4 lead to the theorem 



statement. □ 

6. Reachable configurations on arbitrary graphs 

In this section, we show that the same asymptotic estimate in Theorem |4.8| for 
the number of debt-reachable configurations also holds for the number of reachable 
configurations. Specifically, we have the following theorem. 

Theorem 6.1. The number of configurations reachable from an arbitrary configu- 
ration with c chips asymptotically approaches 



(n-l)\K(G) 

as c approaches infinity. 

We first prove the following lemma. 

Lemma 6.2. If the configuration (c\, C2, C n ) is debt-reachable from the configu- 
ration (c'j, c' 2 , c' n ) then the configuration (ci + di, c 2 + d 2 , c n + d n ) is reachable 
from configuration (c^ + di,c' 2 + d2,...,c' n + d n ). (Recall that di is the degree of 
vertex i ). 

Proof. Alternatively, the theorem statement is equivalent to the statement that if 
two configurations C and C are debt-reachable from each other, then they are 
debt-reachable while going at most di chips into debt on vertex i at any point. We 
shall prove this version of the statement. 

Let x = (x\, x 2l ■ ■ ■ , x n ) be the reduced firing vector for the transition from C to 
C . Consider the following algorithm to reach C from C without going more than 
di chips in debt on vertex i at any point for all vertices i: If there exists an index i 
such that Xi > and such that vertex i is not currently in debt, then fire vertex i 
and decrease Xi by 1. 

We claim that this algorithm terminates only when we reach the configuration 
C (i.e., when all the Xi — 0). To see why this is true, assume to the contrary 
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that vertex i is in debt for all i where x% > 0. Let S the set of indices i where 
Xi > 0. We claim now if that we fire all the vertices in S the appropriate number 
of times (as given by the Xi), at least one of these vertices will remain in debt. To 
prove this, note that if we look at the total number of chips on all the vertices of S 
(which is currently a negative number), this number can only decrease or remain 
the same when a vertex in S is fired. Hence, after firing all the vertices in S the 
appropriate number of times, the total number of chips on vertices of S will be 
negative, and hence there is at least one vertex in S which is still in debt. But this 
sequence of chip- firing moves was supposed to reach the (legal) configuration C", 
which contains no vertex in debt. Therefore we have reached a contradiction, and 
the statement of the theorem holds, as desired. □ 

We can now prove Theorem |6.1| 



Proof of Theorem 6.1 First, note that for sufficiently large c, for any configu- 
ration with c chips, it must be possible to reach some configuration C with at least 
d{ vertices on each vertex i. Assume C has the form (ci, c 2 , . . . , c n ). By Lemma 
|6.2[ the number of configurations reachable from C is greater than or equal to the 
number of configurations debt-reachable from C = (c\ — d\, C2 — cfe, . . . , Cn — d n ). 
This configuration contains c — 2E chips, where E is the number of edges of G 
(since ^2 d% — 2-E) . By Theorem |4.9| the number of configurations debt-reachable 
from C is asympotically equal to (c — 2E) n ~ 1 /((n — 1)\k(G)) which in turn is 
asymptotically equal to c" _1 /((n — 1)!k(G)). Hence, the number of configurations 
reachable from this configuration with c chips is asympotically greater than or equal 
to c™ _1 /((n — 1)!k(G)). However, the number of reachable configurations must also 
be less than or equal to the number of debt-reachable c onfig urations, which is also 
asymptotically equal to c™ -1 /((n— !)!«(£?)) by Theorem 4.9 Therefore the number 



of configurations reachable from this configuration is asymptotically equal to 



(n- 1)\k(G) 

as c approaches infinity, as desired. □ 



7. Towards quasipolynomiality of the number of reachable 

configurations 



In Theorem 4.1 we show that the number of configurations debt- reachable from 
the configuration (c, 0, . . . , 0) is a quasipolynomial in c. It is natural to then ask 
the same question about reachable configurations; does the number of reachable 
configurations exhibit any quasipolynomial behavior? For the case of C3, Theorem 
|3.3| answers this in the affirmative. What about for other graphs Gl 

For C4, computer-aided computation of the number of reachable configurations 
from the configuration (c, 0, . . . , 0) for 1 < c < 100 suggests that this number is a 
quasipolynomial for c > 4. Similar computations for other graphs (C5, Cq, P2, Pa, 
K4) suggests the following conjecture. 

Conjecture 7.1. The number of configurations reachable from (c, 0, 0, . . . , 0) on 

an arbitrary graph G is a quasipolynomial in c of degree n — 1 for sufficiently large 
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While we are unable to prove this conjecture, we prove the following characteri- 
zation theorem for configurations which are debt-reachable but not reachable 

Theorem 7.2. Let C be a configuration which is debt-reachable but not reachable 
from configuration C . Let x = {x\, X2, ■ • ■ , x n ) be the reduced firing vector between 
C and C . Now, let D be a reachable configuration from C with reduced firing vector 
y = (2/1,3/2, • ■ • ,2/n)- Then there exists a unique choice of D with the following two 
properties 

i. For all i, 3/j < Xi 

ii. If (xi — yi) > 0, then vertex i cannot fire in configuration D. 

Proof. Clearly there exists at least one such configuration D, since we can just 
keep on firing vertices v with x v > until there are no longer any legal firing moves 
remaining. Thus, we need only show uniqueness. Assume there are two distinct 
such configurations, D\ and D%. Consider the vertices in D\ that were fired more 
times than the vertices in D2 (without loss of generality these exist; otherwise, swap 
Di and D%). We claim that in configuration D2, we can fire at least one of these 
vertices. 

To prove this, consider the first time in the transition from C to D\ when a 
vertex is fired for the (a + l)-th time when this vertex is only fired a total of a 
times when reaching Z?2- We claim that this vertex is fireable in configuration D^. 
This follows since the only operation that can reduce the number of chips on this 
vertex is firing this vertex, so therefore, the number of chips on this vertex has to 
be at least as large as the number of chips on this vertex when it was fired for the 
(a+ l)-th time when reaching D\ (since all prior moves when reaching D\ have also 
been made at some point in reaching L^)- Then, since this vertex was fireable in 
configuration D±, this vertex must be fireable in configuration D2, as desired. □ 

In essence, this theorem states that every configuration which is debt-reachable 
but not reachable corresponds uniquely to a "most fired" reachable configuration 
and a firing vector entirely supported on the non-fireable vertices of this configura- 
tion. We hope that this characterization theorem could be used to prove that the 
number of configurations that are debt-reachable but not reachable is a quasipoly- 
nomial in c; together with Theorem |4.1[ this would imply Conjecture |7.1| Note 
the similarities between this line of reasoning and the line of reasoning used in 
Theorems 13.21 and 13.31 

In addition to this characterization theorem, we also prove two lemmas exploring 
the relationship between debt-reachability and reachability. The first lemma shows 
that it does not matter what order we fire vertices in when trying to reach a 
reachable configuration. 

Lemma 7.3. Assume configuration C is reachable from configuration C = (ci, C2, . . 
and this transition has reduced firing vector (xi, X2, ■ ■ • , x n ). Then, if Xi > and 
Ci > di, let C" be the configuration obtained by firing vertex i. Then configuration 
C is reachable from C" . 

Proof. Since configuration C is reachable from configuration C, there is a se- 
quence of legal chip-firing moves to get from C to C. Let these moves (in order) 
be firing the vertices v\, Ua, . . . ,v r . Let Vj be the first occurrence of vertex i in 
this sequence (we know that there must be an occurrence of vertex i in this se- 
quence because Xi > in the reduced firing vector). We claim that the sequence 
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Vj, v\, V2, • • ■ , Vj-i, v j+i, ■ ■ ■ ,v r is also a sequence of legal chip-firing moves (and 
since it is a permuation of the previous sequence, it also sends C to C). Note that 
since firing vertex Vj — i can only increase the number of chips on other vertices, 
and since none of vertices v\ through Vj-i are equal to vertex i, the sequence of 
chip-firing moves Vj, Vi, V2, ■ ■ ■ , Vj-i must be a sequence of legal chip-firing moves 
(the first move is legal by the conditions of the theorem statement). Then, because 
Vj, Vi,V2, ■ ■ ■ , Vj_i is a permutation of vi, V2, ■ ■ ■ , Vj, after j operations the configu- 
ration is the same as it was originally after j operations. Since the remaining moves 
Vj+i, ... ,v r have not changed, they therefore still remain legal. □ 

The second lemma (along with its corollary) shows that we never need to fire 
every vertex to reach a configuration. As a consequence, if we can reach a con- 
figuration, we can reach it by firing each vertex the number of times given in the 
reduced firing vector. 

Lemma 7.4. If configuration C is reachable from configuration C — (c 1; c 2 , . . . , c„) 
through a sequence of moves that has firing vector (xi, X2, ■ ■ ■ , x n ) where all Xi > 0, 
then it is also reachable from configuration C through a sequence of moves that has 
firing vector (x\ — 1, X2 — 1, ■ • • , x n — 1). 

Proof. Again, let the sequence of chip-firing moves to get from C to C be firing 
the vertices V\,V2, ...,v r , in that order. Since Xi > for all i, we know that each 
of the vertices 1 through n appears at least once in this sequence. Without loss of 
generality, let the last occurence of vertex i occur after the last occurrence of vertex 
j iff i < j (we can attain this by simply relabeling the vertices). Now, let Vk n be 
the last occurence of the vertex n, and let v^ be the first occurence of the vertex i 
after Vk n for 1 < i < n — 1. For instance, if our sequence was 

1,4,4,2,2,3,1,2,1,1 

then &4 = 3, fci = 7, k 2 = 4, and £3 = 6. We claim that the sequence obtained 
by removing v^ for all 1 < i < n is still a legal sequence, proving the lemma. 

To show this, assume to the contrary that after we perform this operation, vertex 
v p = i is in debt. Let the set S consist of the vertices that occur before vertex v p in 
the original firing sequence that were removed. By the above construction, S must 
include vertex i, since vertex must occur before vertex i. Also, since we remove 
each vertex just once from the sequence, S contains no repeats. 

Now firing vertex j only decreases the number of vertices on vertex i if j = i. 
Since firing all n vertices does not change the number of chips on vertex i, firing 
a subset of the n vertices which contains vertex i cannot increase the number of 
chips on vertex i. Since S is such a set, firing all the vertices in S cannot increase 
the number of chips on vertex v p . But this means that vertex v p must have been 
in debt in the original sequence as well, which is impossible. Therefore our new 
sequence is a legal chip-firing sequence, as desired. □ 

Corollary 7.5. Assume configuration C is reachable from configuration C — 
[c\, C2, ■ ■ ■ , c n ), and let x = [x\,X2, ■ ■ ■ ,x n ) be the reduced firing vector for the 
corresponding transition. Then we can reach configuration C from configuration C 
by firing vertex i exactly Xi times for each 1 < i < n. 
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8. Remarks and open problems 



The most relevant open problem is Conjecture |7.1| The author of this paper 
believes that there should be a proof of this conjecture using Theorem 1 7 . 2 1 and the 



two Lemmas 7.3 and 7.4 Efforts along these lines have not been successful thus 
far, however. 

If Conjecture |7.1| is answered in the affirmative, there are a host of related ques- 
tions that can be asked. For instance, how large must c be for the number of 
reachable configurations to exhibit quasipolynomial behavior? Preliminary inves- 
tigations suggest that c > 2E — n should suffice, but it is unclear if a better 
bound can be found. Similarly, what properties does the period of the resulting 
quasipolynomial satisfy? Note that for C5 it appears that the period of the associ- 
ated quasipolynomial is 10, so the period likely does not have to divide n(G) as in 
Theorem |4.1 1 We can also ask what these quasipolynomials look like in the case of 



cycle graphs C n , extending the result of Theorem 5.4 

Even in analyzing the number of debt-reachable configurations has met with 
more success, there are several questions that remain to be answered. For instance, 
what are the exact periods in Theorem 4.1 (we know only that they divide re(G))? 



For what other graphs (aside from C n ) can we find particularly simple forms for 
the number of debt-reachable configurations? Is there a particularly simple combi- 
natorial interpretation of the numbers bi for small values of i (for instance, 61)? 

Finally, it would be interesting to extend these methods to the case of particular 
classes of infinite graphs, for instance an infinite path or an infinite square lattice. 
Such research could perhaps be relevant to various chip-firing and sandpile model 
problems in statistical physics. 
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